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On the Integration of a System of Differential 
Equations in Kinematics. 

By John Eiesland. 



If we have a rigid system moving about a fixed point and p, q, r be the 
components of rotation, in order to determine the motion completely we must 
integrate the following system of differential equations, 

da 

- W =Pr-yq, 

d(3 



(1) 



which has been treated by Euler, D'Alembert, and more recently by Darboux 
(Theorie des Surfaces, vol. I, Chap. 1-2). Euler expresses the 9 cosines of the 
variable axes of the system in terms of three angles 6, ^i <?> (the so-called 
Euler's angles) which, if known, will determine the motion of the system. 

Darboux reduces the problem of integration to that of a Riccati equation 
by introducing two new variables x and y which remain constant respectively 
on each set of rectilinear generators of the sphere a 2 + /3 2 + y 1 = 1 (Lepons, 
vol. I, p. 22). This method, while very elegant, does not admit of any obvious 
extension to systems with four and in general n variables such as have been 
deduced by Professor Craig and Mr. Hatzidakis,* while the introduction of angles 
corresponding to Euler's angles becomes impracticable for higher dimensions. 

In a note published in Am. Jour. vol. 20, I have sketched an integration 
theory for a system in 4-dimensional space adopting a method employed by Lie 
to the general linear system in two and three variables. The essential feature 

*Am. Jour, of Math., vols. 30 and 23. 
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of this method is a transformation of the co-ordinates by means of the formulae 
Pi _ _ Ps 



«i 



^P? + p! + ---- + pLi + i' 3 ^p! + pl+ •■•• +p 3 »-i + i 

i 



a n 



Vp 2 1 + pl+ ....+ P U+i' 
which is nothing but an introduction of homogeneous variables. The homogene- 
ous system is thereby transformed into a system of generalized Riccati equations 
in n — 1 variables whose coefficients are functions of t as before. 

A remark may here be made concerning this transformation. Consider the 
case of three variables 

ai = v P f + pi + i' a^= vpfTpTTl , a3= V$+JTi' 

If a 1( a 2 , a s be a point on the unit sphere, then, by central projection of the point 
on the plane a 3 = — 1, the corresponding point in this plane will have the 
co-ordinates p x and p 2 , the point of intersection of the a 3 axis with this plane 
being taken as the origin. To a rotation of the sphere there corresponds a 
certain motion of the points on the plane, namely a Non-Euclidean motion whose 
absolute is the circle of + pf + 1 = 0.* 
The transformed system 



-Jr = — 9. + fV + pi (paP — Pi?) -. 



dt 
-&.= p — 9l r + p z (p ai ) _ m ) , 

is, however, not easier to integrate than the original system. The only advan- 
tage we can draw from this transformation is to render the problem more 
intuitive, if we intepret p 1( p 2 , t as coordinates of three-dimensional space, (See 
Lie-Scheffer, Cont. Gr. p. 778), which advantage is lost in the extension of the 
problem to four or more variables. Moreover, the Euclidean Motion on the sphere 
is more intuitive for very good reasons than the motion of the Non-Euclidean 
manifoldness on which the sphere may be projected centrally. We shall there- 
fore not use this transformation, but consider the system in its homogeneous form. 
It is the object of this paper to show how by another method that is capable 
of generalization the integration problem may be reduced to its simplest terms. 

* See Klein, Nicht-Euclidische Geometrie, Vorlesungshefte, I, p. 223. 
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We shall use geometrical language as much as possible and integrate geometri- 
cally in the sense employed by Lie in his Beruhrungstransformationen, and in 
Lie-Scheffer's volume on continuous groups. The analytical treatment may then 
be considered as a verification of theorems that have been made evident from 
geometrical considerations. 

I. 

We shall begin with the case of two variables and consider the equations 

dx x 



dt 
dx 2 



Puffy • 



(2) 

where p 12 is a given function of t. These equations admit of the Euclidean group 
of rotations around the origin of the plane x x x 2 , that is, the transformation 

but since t is a variable we must conceive of the transformation as continually 
changing with the time t. Now since all the circles x\ + x\ = const, are 
invariant curves of the system, the two lines 

X\ — %X 2 , x± — — %x% , 

which pass through the circular points at infinity are common tangents to all 
these circles and may therefore be considered as their envelope ; hence they are 
a pair of integral curves of (2). But any system of the form 

-^- 1 = a 1 x 1 + / tf 1 * 2 , 



dx 2 



' — a 2 *i + &*» 



can be solved by one quadrature if two integral curves are known, hence (2) can 
be solved by one quadrature only. 

This may easily be verified analytically. Introducing the variable -- 1 = a 

x 2 

in (2) the system reduces to the form 
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which is a Riccati equation having a = ±ifor particular solutions and can there- 
fore be solved by a single quadrature ; the general integral is 

x 1 — x sin a , x 2 = x cos a, 

where a=fp n dt. 

We shall now consider the system in three variables 

dx 1 



(3) 



dx 2 , 

—=r. PlS&l r]p23 X 3 ' 

dx 3 

~JT — Pl3 x l P^Sph y 



which is identical with (l), if we put^>i 3 = q, p 1% = r, p n = p. Since the oo 1 
spheres x\ -f x\ 4- x \ = const, are invariant surfaces, if we determine the motion 
of any one of these, say the unit sphere, the complete motion is known ; we may 
therefore limit ourselves to the consideration of this sphere. 

Let it be supposed that we know an integral curve Con the sphere ; we shall 
prove that the integration of (3) may be reduced to that of the system 

w &="*• $=-«*■ $=•• 

Let the initial position of the x u x 2 , ar 3 -axes be such that the « 3 -axis passes 
through a point P on C at the time t . We now transform these axes to a set of 
new but variable axes y^, y 2 , y 3 such that the 2/ 3 -axis describes the curve C; this 
transformation will be a Euclidean one whose 9 direction-cosines are functions of 
t depending on the given integral curve. "We now pass a plane, y 3 = 0, through 
the origin and perpendicular to the y 3 -&xis and choose arbitrarily in the inter- 
section of this plane with the sphere two points Q and B at a quadrant's distance 
from each other. As y 3 moves along C, the two axes OQ = y lt OR = y 2 will 
describe a unit circle in the plane y 3 = 0. But since the transformation is 
Euclidean, the transformed system will have the same form in the variables 
V\i Vz< V3 as the original system (3). Moreover, y 1 = y 2 =. will be a particular 
solution of the system, that is, the system takes the form (4), q. e. d. 

To verify this theorem by analysis we proceed as follows : 

Let the integral curve be given in the form 

(5) Xi = pj (t) x 3 , x 2 = p 2 (t) x 3 . 
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We transform the axis by means of the orthogonal transformation 

Vi = v-vPh ■+■ ^^ + a 13 x 3 , 
(6) - y 2 =za 21 x 1 + a 2S x 2 + a 2 ^c s , 

. Vs = asi^i + a®P z + a^Cs, 

where the a's are the 9 direction-cosines of the new axis. Now since the y 3 -&x\s 
always passes through a point P on the integral curve, y x — y % =. 0, whenever 
x x = p x x 3 and x% = p 2 a; 3 ; these conditions must then be written 

J¥*u + P2ai2 + a ls =0, 
pi«2i + P2«22 + a 33 = , 



_ «31 «32 

Pi ' P2 

a 33 a 33 



or, more simply, 

(7) 

that is to say, a 31 , a 32 , a 33 are known functions of t, viz.: 

(8) a31 = Vp? + p!+i' a32= N/p? + pI + i' a33 = Vpf+ P I+i- 

Introducing now the new coordinates in (3), we find after easy reductions, keep- 
ing always in mind the 6 relations between the 9 cosines a ik 

-fi = {Pi&m —Pi&32 + PzPsi + <*2ia'ii + ass^K + o^a'si] 2/2 

+ |>i 2 a23 — i>i3a22 + PzP-zi — "Vu — olwfhz — a'ss^is] 2/3- 
_|? = — [^laagj — jJiaOgg + jp^aaj — a 21 a' n — a^a'^ — a za a' ls '] y x 

H" [jPl2 a ]3 i>13 a 12 + i>23 a H a 31 a 21 a 32 a 22 a 33 a 23j Vz> 

-£ S = — IPvfhs — Pv&w + Pzfhi — a'sian — a's2 a i2 — as^is] 2/i 

— [i^is — Pi&k + Px&u — aVzi — a '32«22 — a^as] 2/2 . 

where the a' tt are derivatives of the a's with respect to t. This system may be 
written 



(9) 



-J* P 122/2 -P132/3 > 

-jj — — -P122/1 + -*232/3 1 



%8_ 



^ 



= ^132/1—^232/2. 
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which is of the same form as the original system. But, P 13 = P 83 =0; in fact, 
since by (5) and (7) we have 

(10) -^L = -^l = p -^ = ^*- = p , 

it follows that x x = a 31 , <c 2 = a 32 , a: 3 = a^ are a set of particular solutions of (3) 
and we must therefore also have 

<*31 = i>12«32— B3«33i 
<*32 = — VvP-K + P%&33> 
<*33 = #13031 #23032- 

Substituting these values of a' m , cc' 32 and a'^ in P m and P 23 , we find that they 
vanish identically, and our system reduces to the form 

W % = '*.. f = -P^. ^ = o,j...i 

If now we determine the six cosines a lu a ]2 , a ]3 ; a 21 , a 22 ,a 23 in terms of pjandp^ 
we know the function P n , so that this last system may be integrated by one 
quadrature. Having found y v y 2 and y 3 we substitute in (6) and solve for 
Xi, x 2 and a% which will be the general integrals of (3). Since the knowledge of 
an integral curve C on the unit-sphere by (10) amounts to having a set of 
particular solutions given, we may state the result thus : 

If a set of particular solutions o/(3) is known, the system may be integrated by 
one quadrature. This is the form in which Darboux states this theorem (Lepons, 
vol. I, p. 28). 

It only remains to show how to calculate the six cosines a w , a 2< (*= 1, 2, 3) 
in terms of o Y and p 2 , or, what is the same thing, in terms of a 3] , a 32l a 33 , which 
are known functions of t. 

The integration problem such as Buler conceived it consists in determining 
the 9 cosines of the new or variable axes (which are different from the axes 
2/n 3/2> Vz chosen above) as functions of three angles <|>, 6, ^ (see Darboux, vol. I, 
p. 3). The formulae thus obtained being rather unsymmetrical, he adopted 
another system of parameters, introducing the direction-cosines a, (3,v of an 
axis of rotation a and a rotation 10 about the axis. This representation is simpli- 
fied if we introduce homogeneous parameters putting 

A = sin — cos a , B = sin 4r cos /? , G = sin — cos y, Z) = cos-— , 
2 2 & £ 



. »(*£=) 
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so that J 2 + B 2 -f- G 2 + D 2 = 1 . The nine cosines are then as follows: * 

D 2 + A 2 — B 2 - G 2 , 2{AB—CD) , 2(AC~BD), 
2(AB+GD) , D 2 — A 2 + B 2 —G 2 , 2{BC—AD); 
2(AC — BD) , 2(BC+AD) , D 2 — A 2 — B 2 + C\ 

If we put D — ^, ^ = ^, 5 = £ - C= ~X- where A = 1 + a, 2 + ^ + r 2 we 

obtain the parametric representation of Cayley for n — 3.f The table of cosines 
is then 

1+ "k 2 — ft 2 — v 2 v + V o b» — C- 

2 — v + ty l + fi 2 — v* — tf 

A ' A 

v li + Zv o Vp — Z 1 + v 2 — 7? — (i 2 

A 'A ' A " 

Now since these cosines are connected by two relations only, viz.: 

°%i — _ ( <^ _ p 

a 33 ' a 33 

one of these three parameters is perfectly arbitrary and may be put equal to 
zero ; we have then the following table of cosines, putting Jl = 0, 

1 — fi 2 — v 2 2v —2y. 

A A ' A 

— 2v 1 + fi 2 — v 2 2fiv 
A A ' A 

2(i 2{iv 1 + v 2 — (i 2 

A ' A ' A 

where A = 1 + p' + v 2 . We have now by (7) and (8) 

2fi _ _ pi 2fiv _ _ pa l+v 2 — fi 2 

A -a 31 _ Vi+pf + pi , A -as,- Vl+pf + pl' A 

_ 1 

Solving these equations for n and v we get 



v = _PjL l-Vi + pf + pl 

Pi ' ' Pi 



*See Encyclopaedia der Math. Wiss., Bd. IV, p. 204. 
f Pascal, Die Determinanten p. 163. 
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Remark. This choice of % amounts geometrically to choosing the position 
of the moving axes y 1 and y 2 in the plane y 3 = in such a way that the direction- 
cosines of y 2 become a 21 , a 33 , a®, and those of y lt an, — a 21 , — a 31 . 

Expressing the above cosines in terms of p x and p 2 we obtain the table 



P f— pjvi + pf + pj Pl p 3 (i + vi + p! + pi) -pi 

(pj + pSKi +pi+pg * ( P f + pl)v i + pl + p l' vi + pf + pi' 
- PlPa (i + vi+pf + pl) — pi + pXi + pf + pl +p 3 

(p! + pI)Vi + pi + pl ' (pf+pl) Vi+pf + pl ' vi+p? + pi' 

Pl P2 1 

Vi + pl + pl ' Vi+pf + pl ' Vl+pf + p|' 

or, in terms of <x 31 , a 32 and a^ which is a set of known particular solutions, 

<*33 4" a 31 1 a 31 a 32 



1 a.33 1 G&33 

1 2 
a 31 a 32 1 a 32 a 33 



1 a m . 1 033 

a 31 > a 32 



a 31> 
a 32> 



Substituting these values in P 12 of (4) and integrating we find 

j/ 1 = 3csin6>, ^ = jccosm, y s = G. 

where a=./P 12 dt. Substituting the values of y x , y 2 , y 3 in (6) and solving for 
the x's we obtain the general integral 



(11) 



Xi = a33 + <-l x sin a _ JW?2_ xC0Sa + Ca 3X , 

1 a 33 1 0t33 

X2 = ^^_ xs[ua + 1 — <4» - (X33 x cog u + ^ j 
1 a33 1 — a33 



. a 3 = — a m x sin co -j- a^x cos a + ^o^ . 

Suppose now that two integral curves are known. We shall show that the 
system may be reduced to the form 

(12) % = f? = 0, ^- 3 =0. 

7 at at at 

Let the integral curves be G x and C!;. We start from an initial position, 
placing the extremity of the y s -&xis at a point P on G 1 ; on G 2 we take a point $ 
on the y 8 -axis at a quadrant's distance from P. As t varies, P moves along G x 
while Q moves along G 2 , these curves being integral curves. The axis y x will 
describe an integral curve C 3 which is uniquely determined by the motion, when- 



(14) 
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ever C 1 and C z are known. The transformed system must now admit y x = y 2 = ° 
and y 1 = y 3 = as particular solutions, that is, it must reduce to the form (12), 
q.e.d. 

It follows from this theorem that the system (3) may be integrated without 
quadratures, whenever two integral curves are known . 

Analytically we proceed as follows : 

Let the integral curves be 

(13) «i = pi^3. a-2 = P2 a; 3> x 1 =a 1 x s , x 2 == c 2 x 3 . 

The transformation (6) must now satisfy the relations 

' Pi«n + P2«i3 + a 13 = , 

Pl«21 +P2a22 + «23=0, 

ffiaii + tf2ai2 + ai3 = 0, 

. Ol^l + 0'2 a 33 + «33 = ° ' 

Since y x = y 2 = whenever a^ = pj a; 3 and a 2 = p 3 x 3 and also «/j = y 3 = for 
Xj = cr^ and x z = o%x 3 . These four conditions reduce to three; in fact, the 
distance between P and Q must equal V2 for all values of t, the initial position 
having been chosen as explained above. An easy calculation will show that 
we must have 
(15) p^ + p^+lEEO. 

The conditions (8) may be written 

n a 31 . a 32 . - — a 21 n -- a 82 

Pi , p2 — i "l — i °2 — ) 

«33 a 33 a 23 a 23 

which show that a 31 , a 33 , a 33 ; a 21 , a 22 , a 23 are two sets of particular solutions. 
The remaining three cosines are now uniquely determined, so that the table of 
cosines will be 

^2 — p2 _ Pi — a i __ 

Vpi + ff+l *M + 4 + 1 ' Vpf + pi + 1 V7f+ el 4- 1 ' 

p2< T i — Pi^a 



^Pl + pi + 1 ^1 + <^I + 1 
1 



Vp? + pl+l ' VoJ + of+1' V&i + ol+l' 

Pi P2 1 

Vp? + pi + i ' Vpf + pi+i' Vpf + pi + i' 

(That these quantities are direction cosines may be verified by aid of (15) ). 
4 
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Transforming now to new axes y lt y 2 , y s and reducing we find that P 18 , P 23 
and P n all vanish so that the system reduces to the form (12). The general 
integral is now 



X x = 



Ps 



x 1 + 



OS 



*2 + 



Pi 



«8: 



Vp? + o\ + 1 Vcr? + a\ + 1 '^ s/ol + al+1 2 T Vp? + pi + 1 ' 



Pi 



*1 + 



*2 + 



p2 



*3 



Vpf + pi +1 Va? + at +1 l Vol + <r§ + 1 8 Vp| + pi + 1 ' 



,x 3 = 



Pa^i — p i^a 



x x + 



*2 + 



Vp? + ? \ + 1 Vcxf + al + 1 * ' Vof + o| + 1 ' ' Vpf + pi + 1 



x 3 



II. 

We shall now take up the study of the system 



(1) 



dt 
dx, 



= Pl2%2 — Pl3^3 + Pl&i . 
= — JPU&1 + 2>23«3 + 2fc»4 - 



3 = Mi-W+Wi 






which has been obtained by Professor Craig.* Since the hypersphere 
x\ + x\ -f- x\ + »I = const, remains invariant during the motion we may, as 
before, limit ourselves to the unit sphere. The problem of integration is then : 
given the six components of rotation p ik around the six planes of the hexahedron 
of reference, to determine the motion completely. Let a curve C x on the hyper- 
sphere be known, and let the initial position of the axis at time ^ be such that 
the a^-axis pierces the sphere at a point P on this curve. Transforming to anew 
system of variable axes y x , y%,y 3 , y± in such a way as to make P describe the 
curve G h as t varies, the motion of the other three axes will take place in a three- 
dimensional space obtained by cutting the hypersphere by a space y t = passing 
through the origin ; (this space will necessarily vary in its position with t. It 
follows then that after the transformation y x = y t = y 3 = must be a particular 

* " On the displacement depending on one, two and three variables in space of 4 dimensions," Am. Jour, 
of Math., vol. 20. 
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solution of the transformed system. But since a Euclidean transformation does 
not alter the form of the system (1) we must have 

dyi 



(2) 



dt 



— -^12^2 Pvzy$ ' 



—T. P\%H\ + PzsHJz ) 

—jj- = -fWl PisR/z , 

^y±— o 

dt 
Hence the 

Theorem. If an integral curve of the system (l)is knoicn the integration of the 
system may be reduced to that of a similar system in three variables. 

If an integral curve C' 2 of (2) is known, it is evident that a corresponding 
integral curve of (1) is known and may be obtained by transforming C 2 by the 
inverse of the transformation which changes the system (1) into (2). Let this 
transformed curve be C 2 and let us suppose it different from 6^; then two integral 
curves of (1) are known. But we have proved that if an integral curve of a 
system in three variables is known, it may be integrated by one quadrature ; it 
follows therefore that the original system can be integrated by one quadrature, 
whenever two different integral curves are known . 

Since the knowledge of two integral curves amounts to having given two 
sets of particular solutions, the last theorem may be stated thus : 

// two sets of particular solutions of the system (1) are known, it may be inte- 
grated by one quadrature only. 

Since the integration of (2) by Darboux's method* may be reduced to that 
of an ordinary Riccati equation, the above theorem may also be stated thus : 
If one set of particular solutions is knoicn the integration of (l) may be reduced to 
that of a Riccati equation. 

The analytical treatment of the integration problem is similar to that 
employed in the case of ordinary space. 

Let the integral curves be 
(3) x l = Xj {t) x i} x 2 = X 2 (t) x it x 3 = X z (t) x t , 

where the /l's are known functions of t. We introduce new axes y x , y 2 , y 3 , ?/ 4 by 

* Darboux, Theorie dea Surfaces, vol. I, p. 22. 
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(4) 



means of the Euclidean transformation 

' V\ — an^i + cc 12 x z + a ls x s + a u Xi , 

V% = aai^i + a 22»2 + <ha x 3 + <*u x i > 

Vz = a u x i + <*-3z x z + a^ + a34«4 » 

Vi — a«*i + a*^ + a^ + a^ , 
where the a's are the 16 direction-cosines of the new axis. Now in order that 
the point P shall during the motion describe the integral curve (3), we must 
have 3/1 = 2/2 = 2/3=0 whenever x r = /^ x v x % = "k % x±, x 3 = X 3 x if or 

a ll™l "1" a 12™2 ' a 13^-3 "r <*14 == " , 
a 2 i^i + 0122^2 + a23^3 + «24 = , 
<*31^1 + <*32^2 + a33^3 + OC34 = , 

which may be replaced by the following simple ones : 



(6) 



(5*) 



a, 1 = 






a 42 



Aq — 



_ a« 



a« 



that is to say, a-j = a 41 , x 2 = a 42 , ir 3 = a 43 , x t = % is a set of particular solutions 
of the system (1). Introducing the new coordinates into the system taking into 
account that the a's are a system of 16 direction-cosines of which the four 
a «i a 42> a 43> a u are particular solutions, we obtain after easy reductions 



(6) 



d V\ — 



alt 
d}h 



— P 122/2 P 13I/3 ' 



dt 



*•= 0. 

dt 
where P ]2 , P 13 , P 23 have the following values: 

P J2 = (anOgg — a 2] a ffi ) ^^ + (a 21 a 13 — ccua^) £> 13 -f (a^^ — a 21 a 14 ) p u 

+ (ai2«23 — a22«13) P23 + (<*12«24 «8«n) 2>24 + (a24<*14 ~ G^u) 2%4 

"I Otll^l T Cj 2 a 22 -j- 0.1^23 T &14 a 24» 

-P13 = (ai2«3i — an<%) i>i8 + (an«33 — <*3iai3) P13 + (<%ai4 — anO^) Pu 

+ («32ai3 — <*12«33) i>23 + (ai4<*32 ~ V<l&-3i) Pa* + (a33«14 ai3«34) _?>34 

+ a 31 a u + a^^ + a^a^ 4* a 32 a 12 , 

-fW = (a21«32 — a31«22) PlZ + (a31<*23 ~ a33«2l) Pl3 + («34«21 <*>9l<X>u) Pli 

+ (a33<*22 — a32«23) P-x + (a M a 8 — a^o^) ^34 + (a34«23 — o-3s^2i) Pm 

+ a 31 a 21 + 0^22 + <X-33 a 23 + a 34 a 24 • 



(7) 
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That the coefficients P u , P 24 , P u vanish identically follows from the fact that 
a 41 , a 42 , (X43, a 44 is a set of particular solutions, so that we have 

ail = PrAiz — P\&® + Pn&u . 
a« = — i?12«41 + PaPu + Pu a u • 

«« = i>13«« — PiS^ + 2>34<*44 . 
«« = ~i>M«« i>24«42 — ^34«43 ! 

Substituting these values of a' 41) a' 42( a' 4S , a' 44 in P 14 , P 24 and P 34 they vanish 
identically. 

Since the a's may be expressed as functions of 6 parameters and since by 
(5) or (5 1 ) we have determined 4 of the 16 cosines as functions of \, /t 2 and 2, 3 , 
there remain 3 arbitrary parameters which may be chosen entirely at pleasure. 
In the case of a system in three variables we showed how to actually calculate 
the simplest system of cosines in terms of the known functions g-i and p 2 . The 
corresponding problem for a system in four variables offers no difficulties. If we 
express the 16 cosines a (k in terms of 6 parameters s lt s 2 , s 8 , s 4 , s 6 , s 6 and then 
put s ) = s 6 = s ( =Owe obtain the following system of cosines :* 



(8) 



1 — s'j — si — 4 2«j — 2s 2 2« ; 



3 



A A A A 

- 2*1 1 + 4 + 4— _ s l 2^ — 2^3 
"A ' A '" " ' A ' A 

2s 2 2S& 1 + s\ + s\ — s\ 2s 2 s 3 

A A A A 

- 2s 3 — 2s£ 3 2s 2 s s 1 + s\ + s\ — 4 



A A A A 

where A = 1 +sf •+- s| + s\. Now since a 41 , a 42 , a 43 , a 44 is a set of particular 
solutions we have 

2s 8 — 2s x s z 2s 2 s 3 1 + sf -|- s§ — s% 

— a 41 , ^ a® , — ^ a 43 , ^ a u ; 

Solving these equations for s lt s 2 and s 3 we find 

„ — a 42 „ a 43 „ — a u — 1 . 

&1 , *2 > s 3 » 

a ti a « «. 41 

*See Pascal, Die Determinanten, pp. 160-162. 
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the table (8) may now be written 



g|i + «44 — i 
1— a u 

a 42 a 41 



1 — a u 

«43 a 41 

1 — a 44 



a 42 a 41 

1 — a u 

1 — <4i — «44 

1 — a 44 

a 42 a 43 



— a u 

«42 



«43«41 

1 — a tt 

a 42 a 43 

1 — a 44 

1 «43 «44 

1— «44 

<*43 



— a 41 . 



«42 



and we thus arrive at the following orthogonal transformation : 



(9) 



m = < + -u-l Xi + 



y 2 = 



2/4 = 



1 — an 

a 41 a 42 

1 — a 44 
«43 a 4i 



a 42 a 41 

I— a 44 



Xo 



3*1 + 



<4- 



■a 44 



x-. 



1 — a 44 

a 42 a 43 „ 

— " 2 

1 — a 44 " 
+ a 42 a 2 





+ 




«43«41 

1 — a 44 


<c 3 


iC 3 


— 




a 42 a 43 

1 — a u 


3% 




4- 


1 


— «43 


a 44 



a 43 > 

a 44 , 

— a 41 » 4 , 

"T a 42 a3 4 > 

a*3 + oc 43 a? 4 , 



a 44 
+ <x 43 a% + dit^ , 



1 — a 44 
a 41 Xj 

&?/ means of which the system (1) is transformed into the system (6). 

Now suppose we know a set of particular solutions of (6), or, to use 
geometrical language, an integral curve on the sphere y\ + y% + yi= 1, y t = 0. 
Since the first three equations (6) are independent of the last, we consider these 
equations as defining a rotation of the sphere y\ + y\ + y% = const, around a 
fixed point 0. This case has been treated in the first part of this paper, where 
we found that such a system is integrable by one quadrature, a set of particular 
solutions being known. 

Now I say that if a set of particular solutions of (6) is known, a set of partic- 
ular solutions of (1) is known and conversely. The proof is immediate, for, let 
p* 31 , /? 33 , /? 33 0, be such a set; substituting these for y lt y % , y 3 , y± respectively in 
(9) and solving for the oj's we obtain 

„o a 41 T" &u 1 



(10) 



x\ — 



4 = 



<4 = 



x\ = 



1 — a u 




a 42 a 41 

1 — a u 


An 


a 43 a 41 


An 



&1- 



a 41 a 42 
1 — <*44 



Pm - 



^aPm. 



+ 



+ 



1— c& 



a u 



1 — a tt 

a 43<*41 a 

Pa 



P* 



a 43 a 41 

1 — a 44 

a 42 a 43 



P33> 
Ps3> 



1 — a 44 

a^3. 



+ 

+ 



1 — a u 

1 — a ^ a 44 n 
i — a 44 

a 43p33 > 
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which are easily seen to satisfy the relation arf -f xf + xf + xf = 1. Morever, 
since y± = is a particular solution of the last of equations (6) we have from 
the last of equations (9) 

2/4 = = a a x\ + a^4 + a-isA + ^u x l . 

which means that x\, x%, x\, x\ are four direction-cosines on the hypersphere. 
Conversely, suppose given two sets of particular solutions of (1), viz.: 

a 41 a 42 a 43 a 44 J *1 > • C 2 > &8 > •"* » 

and if we transform (1) by means of (9) into the system (6), this system will have 
for particular solutions (3 n , /? 32 , p^, obtained by substituting x\, x%, x%, x\ for 
the x's in (9), that is, we will have 

Vi = An , 2/ a = Ab , 2/3 = ^33, 2/4 = 0. 2- e -d- 

We are now ready to complete the problem. Transforming (6) to a system of 
new axes z lt z 2 . z$ by means of the transformation 

\ = §n±fo-l yx + -&£f-y t -ft,*, 

1 P33 x P33 

„ — P31P32 ,, 1 1 P33 P32 7/ I a .. 

% — ^ o— Vl + 5 5 2^2 T P322/3 > 

1 — P33 A Paa 

% = Pv&l + Ab?S + ftajfo , 

Z 4 = JA ' 

we obtain a new system of the form 

<»> -?r = ? ^ £ = -*- ^ = ' & = •■ 

which may be integrated by one quadrature. The general integral of (1) may 
now be easily obtained. The general integral of (6) is, (compare (11) part I) 

yi _ fl« + P*l —± x sin a - -p&f- x cos a + Cft, , 

1 P33 1 P33 

y s = 1*£* x S in a + * -/*-&** x cos 5 + Cft,, 

1 — P33 * — " P33 

2/3 = — Pa x sin o + ftj x cos w + Cftg , 
?/ 4 = const. 

where w =/P n dt. Substituting these values in (9) and solving for the jc's we 



(11) 
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obtain the general integral 

a 41 a 42p3lP32 



U= Hi + "44-1 As r /%-! 
* ( l-a u ■ l-/?33 (i-a^i-^) 



1 — a« I 



_ Ufl + «44 1 ffglfe _|_ <*41<*42 1 — ^33 — ^32 

1 1 — a u '1—1^33 1 — ««' 1-/^33 

_0430A) CQS - 

1 — a 44 j 



x 2 



+ c»i + const. a 41 . 

<*44 i — P 



— I <*42<*4i &3 + /?3i— 1 i 1— afa — a^ p^A 
U-<V l-# "*" 1 — « 44 "1 — ^ 



1 — an ) 



_ f <*42«4i _ igA _ 1 — af 2 — «« 1 — ^33 — ^1 
(1 — a^' l—^ 



1— «44 



(13) •{ 



1-&3 

+ a *2 a A } 3 cog - 

1 — a^J 



+ cs4 + const. a 42 

X — -1 CX ^ a » /^33 + /3Il — 1 __ «42«43 J#3lAi2 



1 — a u 



i-A 



33 



1 — <*«' 1—^33 

_ l-«!3-«44 /?3i L sin - 
J «44 ) 



"-43^41 _ PsiPsi I « 42<*43 * Al3i^32 

1 — a 44 1 — /?gg 1 — a 44 * 1 — (3 & 



+ car" + const. a 43 



— a 43 a 44 

1 — a 44 



fl 



32 



X COS 0) 



x 4 



=_j a „.&±i 



A + ffs-l a«&«, 



32 



+ a43 f 



:P*31 [ * 



sin a) 



I f a4lP3lP32 «42P3l P32 _i_ „ a ) „ „. Q - 

+ {r-^33 1=^ + ««P32 J X COS 0) 

+ cx 4 + const. a 44 . 

In these formulae the values of (3 31 , ($&, /? 33 are to be replaced by their values in 
terms of a;?, xl, xl, x\ obtained from equations (10). The determination of the 
angle a is the only quadrature involved, otherwise the work has been purely 
algebraic* 

*It should be noticed that the system of Integrals (13) contain 4 constants of integration: k, O, const, 
and a fourth is obtained in the quadrature fP l2 dt. 
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It is now comparatively easy to obtain an orthogonal transformation which 
shall transform (l) directly into (12). In fact, combining (9) and (11) we obtain 
after somewhat lengthy reductions : 

g _ [ 033 + fill— ! _ «44 + all— 1 _ 031032 _ a4i<*42 _[_ /^lgggji ] x 

1 1 033 " 1 «44 1 033 ' 1 — a« 1 «44 J 

I I 033 + 031 — 1 _ «42tt41 i 031032 1 ~ «44 ~ «44 _^_ 031«42«43 j 3. 

I 1—033 '1 <*44 1 033* 1 <*44 1 «44 J 

+ j 033 + 01l — 1 ^ a 4 3a41 _ 031032 _ a42«43 031 (1~ «48 «tt) ] g 

\ 1 033 '1 <*44 1 033 ' 1 — a 4 4 1 — GC44 j 3 



(14) 



033 + 031 1 „ 0310; 



a « 1 TT~ • a 42 "T" 031 a 43 f ^4 



031^43 J 



1-/8. -" l-0 3 

_ j 031032 _ <*44 + «41 1 _|_ 1 — 033 — 032 _ «41«42 1 032<*43<*41 j x 

(1 — 038" 1~ OC44 1—033 '1 0-44 1 ~ «44 J 

{ 031032 a 42 a 4 1 1 033 032 1 a 42 a 44 1 032 a 43 a 42 I x 
1 033 ' 1 — «44 1 033 ' 1 <*44 1 «44 i 

031032 a43«41 I 1 — 033 — 032 a-42«43 _ 032 ( X «43 O-ii) [ 



X-y 



\l (3 m 'l CC44 1 (3s3 'l a 44 1— «44 ) 

1 I r*3\r^32 1 r^33 " /-^32 i Q \ 

T \ 1 q- • «-41 I ; 7> • «42 T P32 a 43 f X i • 

U 033 1—033 J 

Z3 — XjXj "T" aigiBg T "^3"^3 I" ^4834 ) 

2 4 = a^ + a i2 x 2 + a^ + 0144*4 • 

In the first two equations /? 31 , Z^, 033 have the same values as in (13). This 
transformation exhibits a system of direction-cosines depending on the 8 para- 
meters a ilt a i2 , oc.43, 0144, x\, z%, a^, x\ between which 3 relations exist. We have 
thus incidentally solved the problem : of a set of 16 direction-cosines 8 are given ; 
it is required to calculate the remaining 8 in terms of the given parameters. The 
solution of this problem is equivalent to reducing the integration of (1) to a single 
quadrature, when two sets of particular solutions are known. It should be noticed 
that this problem can be solved in an infinite number of ways, the simplest set of 
cosines being obtained by puttings = s B = s 6 = 0. There exist then a> 3 Euclidean 
transformations which will effect the reduction of the system (1) to the form 

dzi _ P _ <feg _ p _ dz 3 —n dZa — n 

-dJ- F ^> -df- — ^ 1 ' -w- ' Tt u - 

The case where three sets of particular solutions are known (or three integral 
curves) requires no special consideration ; a Euclidean transformation may then 
5 
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dt 



0. 



be found which will transform (1) into the system 

fyi _o #2—o d &» — o 
dt ~ * dt ~ ' eft ' 

which requires no quadrature. If a a , a®, a^, a u ; cc 31 , a&, a^, a^; a 21 , a^, a^, a 24 
be the given sets, the transformation may be written down at once as follows : 

V\ — Al*l + A*«2 + #13*3 + Al«4 » 

2/2 = asa^i + <*-w»h + a^ + a^x* , 

Vz — affl^i -r a^Xg + 0133X3 + 0^X4 , 

Pi — a 4 iXi + a^xjs + a^ + a M x 4 , 
where D n , Z> 12 , Z> 13 , D u are the minors of a n , a 12 , a^, a u in the determinant of 
the 1 6 cosines I a n , a*®, a m , a u | . 

III. 

There remains yet the problem of generalizing the preceding theory of 
integration. Mr. Hatzidakis has obtained the general system denning the rota- 
tion of a rigid system in n-dimensional space around a fixed point.* This system 
may with slight modifications affecting only the signs of the coefficients p ik be 
be written f 

dx x 



(1) 



dt 

2 _, - 



dt 



— VvF<2, + Pl&3 + + PuPn . 

P1&1 + Pis&z + + Pu^n » 



dx n 
dt 



= —PlnXl—Pzn^Z — 



'Pn — ln^n — l' 



It will be easy to prove as in the case of four dimensions that the integration 
of this system may be reduced to that of a similar system, but of one dimension less, 
whenever a set of particular solutions is hnown. The given integral curve being 
G x we transform to a new set of axes yi,y 2 , .... y n in such a way that the y n -&xis 
pierces the hypersphere 2x| = 1 at a point P on the given integral curve. In 

* Am. Jour, of Math., vol. 23. 

tSee my article "On Nullsystems in Space of 5 Dimensions and their Relation to Ordinary Space," Am. 
Jour, of Math., vol. 25, p. 105. Them's in this article are constants, while here they are given functions of t ; 
moreover, n is not restricted to an odd integer; n may be even or odd. 
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the transformed system y 1 = y^-=.y z z=. . . . . y n _ 1 = must be a particular solu- 
tion, so that the last equation reduces to the form -J^ = 0, and the first n — 1 

equations will define a rotation of the sphere y\ 4- y% + . . . . y\- X = C, y n = 
about the fixed point 0. In general, if m integral curves are known, we may 
transform the system by means of a Euclidean transformation in such a way that 
y n describes G ly y n _i describes (7 2 . . . . and y n _ TO describes G m . . . . The extrem- 
ities of the remaining axis will move along paths on the n — rn dimensional 
sphere y\ + y\-\- yl-m— 1» yn-» + i = #.-» + a = y n = 0, and the sys- 
tem (1) will be reduced to the form 

"J 1 = Pl2Vz + Pim + + Pln-rnVn-m, 

t, = PlzVl + -^23^3 "1" + P^n-rnVn-my 



(2) 



ay n -m p p p 

jj J In — mi/1 ■*■ 2n — mif2 - 1 n-m-l.n-m^n-m-li 

= o, 



dt 

dy n - m +i 



dt 



^l*. = 0. 

dt 

In particular, if n — 2 integral curves are known the system may be reduced to 
the form 

%i __ p , dy 2 __ p dy 3 _ Q dy n _ 

and may therefore be integrated by one quadrature. 

We shall show how this reduction may be effected by means of a series of 
Euclidean transformations. The problem consists in determining a set of 
n 2 cosines in terms of a given set of particular solutions, n being any given 
integer > 4. 

Let the integral curve be 

^1 = fl'EiH ®Z = '- Pz^nt > &n — 1 = pn-l^n* 
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Since on the unit hypersphere x\ + x\ + x\ = 1, this is equivalent to the 

knowledge of a set of particular solutions 

pi p% ?n-l 1 

VT+2p„ a _7 Vl+2p„ 3 _x' Vl + SpJ-V Vl+Sp^/ 

We employ the orthogonal transformation 



(3) 



" ?/! = «!!»! + a 13 x 2 + + a ln x n , 

y 2 = a a x x + a 33 aJ 2 + + a Zn x n) 



2/„ = a nl Xi + a n3 a% + + a w „a„. 



In the same way as before in the case of n = 4, y x = y % = • • • • =y n -i must 

vanish, whenever a^ = p^, a- 2 = p 2 x n , se n -i = p«-i«». and hence the 

following conditions must be satisfied : 

ajipi + (X12P3 + + ai«-ip»j-i + a lTC = 0, 

«2lPl + a 22p2 + H~ a 8»— lP» — 1 + a 2« — ° ' 



«»i-llPl + a « — 12P2 + + tt«-l»-lP»-l + a n-ln 0, 

which » — 1 conditions may be reduced to the form 

(4) pi — — . p2— — > f >n - i ~^r^' 

that is, a Bl , a n8 , a„„ must be the given set of particular solutions. In order 

to obtain a transformation that will reduce the system (1), which we shall denote 
by A w , to a system A„_ x , we must find a table of cosines a ik expressed in terms 
of the n cosines a ni , (i= 1, 2 n). Since these n cosines are determined as 

functions of p 1( p 2 . . • • p„_i, there remain w,w ~ 1 — (n— l) parameters which 

may be disposed of at will. Let the — — -^ — parameters be denoted by 

s 1) s 3 , .... 8..-1. We now put all except the first n — 1 of these equal to 
zero and calculate the » 2 direction-cosines in terms of these parameters by the 
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ordinary method. (Pascal, Die Determinanten, pp. 160-161). We find after 
easy calculations the following table, 

1— gf —4 — ••••gjt-i 2*i_ —2s, (—1)^^ 

A ' A ' A ' A 

2 Sl 1 + 4+ .... +*;_ x — si 2s lSz (_i)— i 2glg ._ 1 



A ' A ' A ' A 

2s, 2^ l+4 + 4 + 4 + ---- + ^-x~4 (- If-Hs^^ 

A ' A ' A ' A 



(- l)"- 1 2 g „_ 1 (- ly-*^,-! (— lj^ggjfa 

A ' A ' A ' 

1 ~r s i T g a" " i ~ • • • • + *«— i — *»— i 
A 

n-l 



where A = 1 + 2 S » • ^ e nave * nen 
i 

X a »l) A n2 ' 



1 + sf + + s- n _ % — g;_! 



solving which for s lt s 2 , .... s n _ 1 vre obtain 



„ a n2 a «3 „ — i a «4 „ a nre 



i °3 — i > °n — 1 



and we have the following transformation : 

yi = oA + a^-l ^ + ^a.^ + .... + ^1°^^ +(_ i) XiXb) 

2/2 — —^ — »ii ; x a + "- fXT -a; »-i t v V a «2*» 

/|-\ J ••• U<nn *■ "-ran x T "'nn 



. ^ n — Otnl^l T a n2 x 2 H" ■ • •' *»n — 1*8 "f" ^jin 33 *! 
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which transforms (1) into the form 

-%f = P nV% + PnVz + + Pm-iVn-i, 



(6) 



dt 
dy 



P\n-\V\ Pzn—ilJi P n— in-lVn — 2> 



dt 



= 0. 



By a method analogous to the one we employed in the case of n = 4 we may 
prove that if a set of particular solutions of (6) is known, a set of solutions of 
(1) is also known and conversely. Treating the system (6) in the same manner 
as (1) we may reduce it to another system the integration of which is equivalent 
to that of a system of dimension n — 2. The Euclidean transformation trans- 
forming A„ into a system A„_j is found by combining the two transformations 
transforming A„ into A„_i and A n _x into A„_ 2 . In general, calling the 

successive systems A n , A„_!, A„_ m , and the corresponding transformations 

$i>#.-i, • • • /S„_ OT + i,wemayby means of a transformation S— S n S n _x . . . S n _ m + 1 
transform A„ into A„_ m . In particular, if m = n — 2 the transformation 
S 3 will transform A„ into the system A 2 or, 



S= S n 8 f 



dzi _ p dz 2 _ 



Pl2 Z l> 



dZ; 



3 



dt 



0, 



dz. 



dt 



? = 0, 



which may be integrated by one quadrature. 

It appears from the preceding investigations that the problem of reduction 
of A„ to A n _ m amounts to this: Given m .n cosines as functions of t, to calculate 
the remaining n? — nm in terms of the given functions. This is a purely algebraic 
problem, involving only the combination of m — 1 known transformations 



~»» "» 



S n . 



where S n is of the form (5). There exist 



j n -\, •■•• o,_»4.i wucie u„ 10 ui me lunii ^t>;. xucic caiisi 00 2 

transformations which will transform A„ into A„ _ ! of which the transformation 
(5) is the simplest, corresponding, as it does, to the case where — TT — ' n — 



of the parameters S x . . . . S»- »-» have been put equal to zero. The transformation 
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S— S n , S n _i S n _ m + 1 contains (m — l)n * — j£ — / determinate para- 



meters which are given functions of t; it follows therefore that in general there 
exist n - w ~ 1 - | n . m - 1 - «M* + «0 j = (* -m)(n-m- 1) arbitrary 

parameters, so that there are oo ( """" )( "~™ -1) transformations which will reduce 
A n to A„_ m , »i sete of particular solutions being known. In particular, if 
m= n — 1 we find, since now (n — m) (n — m — 1) = 0, that the transformation 
is uniquely determined by the n — 1 given sets of solutions, which is indeed 
obvious, since n — 1 sets of direction-cosines determine the remaining set of n 

cosines. In the transformation S all the v* m ) \ n m — ) arbitrary para- 
meters have been put equal to zero which insures that S is the simplest 
transformation possible. 

If we compare the determinant of the system (14), (part II) with the 
determinant of the systems (9) and (11), we observe that the first is the product 
of the other two. In general, if we write the successive determinants of the 
transformations S n , S n _ 1 S n _ m + 1 as follows : 



D n = 



« • • • 


a ln 




„(»-l)„(n-l) 
a ll a 12 


„(»-l) 
a ln-l 





<4« 


«2n 




a 21 a 22 


„(«-l) 
a 2n-l 





• : 


a (n) 

"-Tin 


. A,-i = 








<*<«> • • • 


„(n-l)„(n-l) 
<*n-ll<*n-12 • • • 


„(n-l) 
a n — ln- 


-.0 








... 


. 


1 




ag" 


">+l)„(»— m+1) „<«-m + l) 

a 12 .... aj n _ m + i 


.. 


. 




<~ 


(*22 


„(n -m + 1) 
.... l*2n- m+1 


... 


. 



-t/ n — m + 1 



(n-m + 1) (n-m + 1) 
<*ra-m + ]l <*»» — m + 12 



(n-m + 1) /-, 

• • t *n-m + ln-m + l " • •• 



1 

and from the product D = D n D„_ 1 . . . . D n _ m + 1 we obtain the determinant of 
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the direction-cosines in the transformation S which transforms A„ into A„_ m . 
Let now m sets of direction-cosines be known particular solutions and let these 
be written 



x, 



«i> 






• X n 



•Hi ^n — 12> 



•^n-ln' 



»n- 



-11' X n-m + l% 



X n -m + ln 



In the transformation S„ which transforms A„into A„_j the last row of direction- 
cosines are identical with the first set of given solutions (see (5)), and all 
the remaining cosines in D n are expressible in terms of these solutions. In 
D n -! the n — 1 st row are particular solutions of the system A„_ x and the 
remaining cosines are expressible in terms of these ; but by means of the 
transformation S' 1 we may express a'^Zj 1 ] oSlzH . 



•^n - In i 



a%Z\\ _i, in terms of the 
so that all the elements of 



second set of solutions x° _ u x\ _ 12 . 

D n _ 1 will be expressed in terms of these functions. Continuing the reasoning 

in this way, we find that the elements of D n _ m + 1 may all be expressed in terms 



x n _ 



of the last set of given solutions, viz.: x° n _ m + n , xl_ m+lz . 

follows that the product D will be an orthogonal determinant of the form 



■In- 



It 



D = 



A 



21 



-^12 

A 



22 



L ln 



A, 



2n 



A n -ml A„ 



.A n 



X; 



n-ro + ll a3 n-m + 12 



X, 



n-m + ln 



x n \ 



x, 



nZ 



x: 



where the A's are functions of the as's in the last m rows, 
determinant reduces to the form 



Ifm = 



n 



1 the 



D — 



A\\ A 12 . . . ■ A x 



x. 



21 X W 



X, 



2n 



X n l X ni 



where A n 



A ln are the minors of a u , xf 2 .... xf n in the determinant 



c ll> X 22> 



X n 
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If m = n — 2, we have the determinant 



D = 



A»A 



12 ■ ■ 
A n A 22 . . 
~o „o 

x 31 • c 32 • • 



.A 
.A 



In 



2n 



x n\ x n"S • • • • X n 



and the general integral of (1) may be written down at once. We have 

x x — A n x sin u +• A 21 x cos a + Gjx9 ai + G<p\ x + + r „_ 2 x" 1 , 

x 2 = A 12 x sin w + A^x cos t> -(- Cia-g, + C 2 x\ 2 + + C n _ 2 x° ni , 



x n = A ln x sin 6) + J. 2n » cos co + Cia&, + C&P in + + C n _ 2 x\ 



o 

H7i i 



where a =fP 12 dt, and G lt G 2 . . G„_ 2 are integration constants. A n , A n . . A ln 
and J. 2 i, J^ .... A% n are all calculable in terms of the n — 2 given sets of 
particular solutions by a purely algebraic process as was explained above. For 
m = n — 1 we have the general integrals : 

x i = C 1 A n + (Pxl, + G 3 x° sl + G n x» nl , 

X 2 == 'l-a-12 "I" ^2^22 "I" VB#88 ~T" ^n X n2i 



X n — ^1-^-ln "T ^% x in T ^> n X nn , 

where A n , A n .... A ln are the above-mentioned minors and G 1% G 2 . ■ ■ ■ G n are 
the n integration constants. 

The integration problem is thus reduced to its simplest form. It may be 
stated thus: If m sets of particular solutions of the system A n be known, the 
integration may be reduced to that of a system A n _ m . If this last system be 
transformed by means of the transformation 

Vl _, % y% _, x 2/n-m-l 



Vn-r 



= 6)1, 



Vn 



y»- 



w »-m-l> 



into a generalized Riccati system in n — m — 1 variables, the above statement 
may be put in the form: If m sets of particular solutions of A n be known, the 
system may be reduced to a system of Riccati equations in n — m — 1 variables, 
6 



42 Eiesland : Integration of a System of Differential Equations in Kinematics. 
For m = « — 2 we find the Riccati equation 

£ = «#+-■>. 

which is integrable by one quadrature. 

It should be noticed that the Riccati equation here obtained differs from the 
one obtained by Darboux's method in the particular case n = 3, m = 1 ; this 
equation is (Theorie des Surfaces, vol. I, p. 22) 

da _ . , p 13 — ip i3 , p 13 + ip 23 . 

and has no analogue for n >■ 3. It was the failure of this method when applied 
to the general system that led me to search for a method which would be 
applicable in general, and it has been shown that this method is chiefly based 
upon orthogonal, or more explicitly, Euclidean transformations. 

U. S. Naval Academy, December 27, 11)04. 



